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The stability of z/ = 1/3 Fractional Chern Insulator (FCI) phase is analysed on the example of 
checkerboard lattice undergoing a transition into Lieb lattice. The transition is performed by the 
addition of a second sublattice, whose coupling to the checkerboard sites is controlled by sublattice 
staggered potential. We investigate the influence of these sites on the many body energy gap between 
three lowest energy states and the fourth state. We consider cases with different complex phases 
acquired in hopping and a model with a flattened topologically nontrivial band. We hnd that an 
interaction with the additional sites either open the single-particle gap or enlarge the existing one, 
which translates into similar effect on the many-particle gap. Evidences of FCI phase for a region 
in a parameter space with larger energy gap are shown by looking at momenta of the three-fold 
degenerate ground state, spectral flow, and quasihole excitation spectrum. 


Recent work on Fractional Chern Insulators (FCI) 
as a lattice version of Fractional Quantum Hall Effect 
(FQHE)ii^ without a need of Landau levels has attracted 
significant attention^ri^. Those are many-particle exten¬ 
sion of Chern insulatorsi^ - systems which exhibit integer 
quantum Hall effect without magnetic field and were re¬ 
cently realized experimentall} d^d^ . FCI are particularly 
interesting because they can mimic Landau level physics 
and may provide a more convenient way of conducting 
experiments on FQHE, as they can exist in higher tem¬ 
perature and would not need high magnetic fields^. FCI 
can also depart from Landau level physics, which hap¬ 
pens e. g. for bands with Chern number higher than 
one, where new forms of FCI states can arisei^”— . 

Experimental realizations of FCI phase were proposed 
in different systems including cold atoms^S or molecules 
in optical lattice s^^i^^ , graphene^^Ti^^, arrays of quantum 
wires 2 S, transition-metal oxide heterostructures 2 ii 2 ^, or 
strongly-correlated electrons in layered oxides^^Ti^. 

Initially, it was proposed that FCI should exist on 
topologically nontrivial flat band models^i^ii^i^“— . Sev¬ 
eral lattice models with quasi-flat topologically nontrivial 
bands have been shown numerically to exhibit FCI phase, 
including checkerboard^”—honeycomt£i^, square^, 
triangular—, and Kagome lattices^. Numerical evidence 
for analogs of a number of FQHE states, including Laugh- 
lin I/tt^, CF hierarchy^i^ and non-Abelian Moore-Read 
and Read-Rezayi state a^^i^^i^^ was found. For bands with 
higher Chern number, states with no direct analog in 
FQHE were found, some of which exhibiting non-Abelian 
statisticsi^diiiS. 

To prove existence of FCI in torus geometry for fill¬ 
ing p/q one should show q quasi-degenerate ground 
states^”—, which flow into each other and do not inter¬ 
sect with higher states when one flux quantum is inserted 
through a handle of the torus^”— , and obey the momen¬ 
tum counting rules^i^. These rules need to be satis¬ 
fied also for quasihole excitations^i^. Alternative meth¬ 
ods of proving FCI existence include many-body Chern 


number— and entanglement spectrunt^i^”— . 

There are several criteria which allow to find systems 
which can host FCI phase. First, the flatness ratio (a 
ratio of magnitude of band dispersion to the energy 
gap) needs to be low, to maximize the effect of inter¬ 
action. However, this criterion has proven ambiguous, 
as the single-particle dispersion can stabilize the FCI 
phas 0 ^ d^d^~ i^^ and interactions far exceeding band gap 
do not always lead to destruction of FCI^. Secondly, 
in the limit of long wavelength and uniform Berry cur¬ 
vature, the projected density operator algebra resembles 
the Girvin-MacDonald-Platzman algebraS for a Landau 
level. In consequence an energy band needs to have 
nearly-flat Berry curvature to host FCI phas o^^d^ . Also, 
a third criterion, based on Fubini-Study metric was pro¬ 
posed recentl}^”— . However, clear conditions for FCI 
existence are not perfectly understood. 

In this work, we want to investigate how the stability 
of FCI on a given lattice is affected by introducing an in¬ 
teraction with extra lattice sites. We consider a checker¬ 
board lattice which transforms into a Lieb lattice^”— 
when a second sublattice is introduced into the system, 
controlled by on-site staggered potential. We investigate 
the transition between two lattices in the context of FCI 
phase for spinless particles for 1/3 filling. For finite-size 
systems in a torus geometry, we analyze the influence of 
the interaction between the two sublattices on the many- 
body energy gap between three lowest energy states and 
the fourth state. For a specific choice of parameters cor¬ 
responding to an area of larger energy gap, we search 
for signatures of 1/3 Laughlin-like phase. Three lowest 
energy states (a three-fold ground state manifold) are 
analyzed with respect to (i) their momenta, (ii) the en¬ 
ergy gap to excited states for different systems sizes, (iii) 
spectral flow. Also, the quasihole spectrum and its mo¬ 
mentum counting is investigated. Our results suggest 
the existence of FCI phase with a stability supported by 
the interaction with extra lattice sites. The paper is or¬ 
ganized as follows: in Section I we describe the lattice 
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model, Section II contains a single particle analysis, in 
Section III many-body effects are investigated, and in 
Section IV we conclude our results. 


I. MODEL 




A face centered 2D square lattice called a Lieb lattice 
is considered, shown in Fig. [U^a). The lattice can be 
divided into two sublattices A and B, distinguished in 
Fig. ma) by red and blue colors. We use tight-binding 
Hamiltonian 


H = t cj Cj + A cj Cj + 

i'ij) ((*d>> 

+ Vst — Vst ( 1 ) 


ieA 


ieB 


where in the first term () denotes summation over nearest 
neighbors with the hopping integral t, the second term 
is a next-nearest neighbors term denoted by (()) with 
hopping amplitude A and an accumulated extra complex 
phase (/)ij = ±0 when going clockwise and counterclock¬ 
wise, respectively, and Vgt is a staggered sublattice po¬ 
tential. We note that for 0 = 7r/2 the second term corre¬ 
sponds to Kane-Male spin-orbit couplingS^, and (p = tt /4 
was considered for checkerboard lattice in Refs^^i^. In 
the latter case, extra hoppings were added to open the 
gap and fiatten one of the bands; they are shown as t2 
and ts in Fig. [Hb). Based on Refs^i^, these hoppings 
have values t2 = ^3 = 2+2V2 ‘^ transition be¬ 

tween a Lieb lattice and a checkerboard lattice is driven 
by tuning Vst to infinity. In this case, lattice sites rep¬ 
resented by red color in Fig. Hb) are decoupled from 
sites represented by blue color, and systems consisting 
of sites of different colors can be treated independently, 
with blue color sites forming a checkerboard lattice. A 
systematic analysis of this transition will be presented in 
next Section. 

Many-body effects are studied using density-density in¬ 
teraction of the form 

V = Vnn V niUj + Vnnn V niUj, (2) 

(*d) ((«d)) 


Vt < 2A w = 2A w > 2A 

(b) 

FIG. 1 : (a) Structure of Lieb lattice. Red and blue atoms 
belong to sublattices A and B, respectively. Solid black 
lines denote real nearest-neighbour hoppings, arrows de¬ 
note complex second-neighbour hoppings. Other solid 
lines denote further-neighbour hoppings used to fiatten 
the middle band. t2 hopping connects the second-nearest 
neighbours within B sublattice, if an A atom is between 
them. Otherwise, the hopping is —^2- ts hoppings con¬ 
nect third-nearest neighbours within B sublattice. Grey 
ellipses denote interaction parameters, (b) Topological 
phase transition in Lieb lattice. 

where Nx{Ny) is a number of unit cells in x{y) direction. 
We consider 1/3 filling of the middle band which corre¬ 
sponds to N = Nx' Ny/3 particles in the system. Due to 
a translation symmetry and momentum conservation of 
two particle Coulomb scattering term, many-body eigen¬ 
states can be indexed by total momentum quantum num¬ 
bers Kx and Ky^ which are the sum of the momentum 
quantum numbers of each of the N particles modulo Nx 
and Ny, respectively. 


where rii is a density operator on site i, and Vnn and 
Vnnn are interactions between first and second neigh¬ 
bors, respectively. We will focus on correlation effects 
within the middle band, so the Hilbert space is trun¬ 
cated, containing states from this band only. The lower 
band is considered as completely filled. Also, a fiat-band 
approximation is used neglecting the kinetic energies. 
We note that middle band states are localized mostly 
on one sublattice (indicated by blue color in Fig. IHb)) 
even for low Vst, as long as it is topologically nontriv¬ 
ial. Therefore, the leading term in Eq. m is between 
second-neighbors, Vnnn- All calculations are performed 
for finite size Nx x Ny samples with a torus geometry. 


II. SINGLE PARTICLE ANALYSIS 

The unit cell of Lieb lattice consists of three sites giv¬ 
ing three energy bands after diagonalization of Hamil¬ 
tonian given by Eq. [TJ A band structure in the sim¬ 
plest case when only nearest-neighbor hopping integrals 
t are included has the lower and upper bands touching 
each other in the middle of energy spectrum at energy 
E = 0, where the perfectly fiat third energy band is 
present^. Two dispersive bands are almost equally lo¬ 
calized on both sublattices, while the fiat middle band is 
almost fully localized on a sublattice indicated by a blue 
color in Eig. [D^a). We next introduce the second term 
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from Eq. [T]with 0 = 7r/2. The energy gap opens and the 
lower and upper bands are topologically nontrivial with 
Chern numbers C = — 1 and C = 1, respectively, and the 
middle flat band is topologically trivial with Chern num¬ 
ber C = 0, as shown in Fig. [IJb) on the left. Following 
Zhao et al^^, the topology of the energy bands can be 
changed by introducing a staggered sublattice potential, 
i.e. the two last terms in Hamiltonian given by Eq. [TJ 
Increase of Vst leads to bending of the middle band. At 
a critical value of Fgt = 2A, the middle and lower bands 
touch, the band structure shown in the middle in Fig. 
m- At this point a topological phase transition occurs. 
For Vst > 2A, the lower band becomes topologically triv¬ 
ial with Chern number (7 = 0, while the middle band 
becomes nontrivial with Chern number (7 = — 1, what is 
shown in Fig. [Hb) on the right. Similar transition occurs 
for 0 = 7r/4, but at the value Vst = and at Vst = A 
when t 2 and ts are considered. 

Two energy gaps are indicated in Fig. m on the 
right, Eg^ between two topologically nontrivial bands, 
the upper ((7 = 1) and the middle band ((7 = — 1), 
and Eg^ between topologically nontrivial middle band 
((7 = —1) and topologically trivial lower band (C = 0). 
We investigate a magnitude of these gaps as a function of 
model parameters. In Fig. [2fa), a schematic evolution of 
the energy bands as a function of a staggered sublattice 
potential Vst for A = 0.2 and 0 = 7r/2 is shown. With an 
increase of Vst increases Eg^ separating two topologically 
nontrivial higher energy bands from the lower band. This 
corresponds also to decoupling of a sublattice indicated 
by a red color from a sublattice indicated by a blue color 
in Fig. [D^a). In a limit of Wt ^ oo, two sublattices are 
completely decoupled and the Lieb lattice transforms into 
the checkerboard lattice, (blue sites in Fig. [Hb)). At the 
same time, the energy gap Eg^ between two topologically 
nontrivial bands from the checkerboard lattice decreases 
monotonically to zero. A map of a magnitude of the en¬ 
ergy gap Eg^ as a function of the staggered sublattice po¬ 
tential Vst and A for 0 = 7r/2 is shown in Fig. EJb). The 
staggered sublattice potential Vst is varied from Vst = 0 
to Vst cxo, which can be performed by introducing a 
parameter s given by a formula Vst = 4tan(s7r/2), where 
s changes in a range of values s = (0,1). For an isolated 
checkerboard lattice corresponding to s = 1 (I^t ^ oo), 
Eg^ = 0. An introduction of finite Vst opens the energy 
gap Eg ^. 

For sufficiently high Vst the energy gap is a direct 
gap in M point of the Brillouin zone, with magnitude 
= \/4t2 + Fj - Vst • Below Vst = ^ - 2A (white 
line in figISfb)) the bottom of the highest band is located 
at r point, therefore Egi is an indirect gap of magnitude 
4A. We note that the bandwidth of the middle band in 
the topologically nontrivial region is also 4A, so the flat¬ 
ness ratio of the middle band is < 1. The energy gaps 
for phase 7r/4 show similar behaviour, although closed- 
form expression for Egi for high Vst cannot be obtained. 
On the other hand, for 0 = 7r/2 and Vst oo the bands 
touch at the whole boundary of the Brillouin zone (hence 



FIG. 2: (a), (c) Evolution of band structure of Lieb lat¬ 
tice in function of Vst for fixed A = 0.2. (b),(d) Maps of 
single-particle energy gap Egi depending on A and stag¬ 
gered potential Vst^ parametrized by Vst = Atan(s7r/2). 
The top and bottom rows correspond to 0 = 7r/2, and 
(j) = 7r/4, respectively. The white line in (b) shows the 
border between direct (above the line) and indirect (be¬ 
low the line) gap between the upper and middle bands. 


the energy of highest band at both M and K points 
asymptotically approach the top of the middle band), 
while for 0 = 7r/4 the gap is closed only at M point. 

If additional hoppings are included for 0 = 7r/4 (Figs 
[2jc), (d)), the top of middle band is not located in any 
high-symmetry point, therefore Egi can be obtained only 
numerically. In Fig. [2fc) we show dependence of energy 
gaps on Vst for A = 0.2. Similarily to the previous case, 
Eg 2 increases to infinity with increased Vst- However, 
contrary to the previous case, in the Vst oc limit Egi 
remains finite (as was noted in Re^, additional hoppings 
open the gap for checkerboard model). As shown in Fig. 
Hd), the value of this gap depends on A, which is the 
only single-particle energy scale in Vst cxo limit. We 
note that Egi for given A has maximum for finite Vst 
(Fig. [Ifd)), e.g. at s ^ 0.3 for A = 0.2. Therefore the 
additional atoms increase the energy gap, which may be 
beneficial for stability of FCI states. 


III. MANY-BODY RESULTS 


A. The transition between Lieb and checkerboard 
lattices 

The staggered sublattice potential Vst controls the en¬ 
ergetic distance between sites forming a checkerboard lat- 
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tice (indicated by blue color in Fig. [U^b)) and extra sites 
(indicated by red color in Fig. [D^b)) introduced to create 
the Lieb lattice. Analyzing an existence of a Laughlin- 
like phase during the transition between two lattices, we 
look at the magnitude of the energy gap between 3-fold 
degenerate ground states and the fourth state. We per¬ 
form calculations on (4 x 6) torus for interaction param¬ 
eters Vnn = 1-5 and Vnnn = 1- Despite the flatness 
ratio not exceeding one, flat-band approximation^*^ is 
applied as a first approximation, to focus only on the 
effects of interaction, neglecting effects of single-particle 
dispersion and mixing with other bands. In our calcula¬ 
tions, we assume the lower band is completely filled and 
the middle band is filled in 1/3. We have verified the 
validity of neglecting of excitations from the lower band 
checking that they do not significantly affect the many- 
body energy of the three lowest states. We only noticed 
some effect of electrons from the lower band close to a 
single particle topological phase transition, where results 
should be treated tentatively. 

We first consider the situation for 0 = 7r/2 in the sec¬ 
ond term of Hamiltonian given by Eq. [H Fig. [3] on the 
left shows a map of the energy gap as a function of A 
and a staggered sublattice potential, represented by the 
parameter s. A single particle topological phase tran¬ 
sition is marked by a white line in the graph, with a 
topologically nontrivial region above the line. Opening 
of the energy gap coincides with single-particle topolog¬ 
ical phase transition for Vst = 2A, similarly to results 
from Refj^. Within a topologically nontrivial region the 
energy spread S of 3-fold degenerate ground state does 
not exceed S = 0.015. Therefore, in a major part of this 
region 3-fold degenerate ground state separated by the 
gap is clearly seen in the energy spectrum. Values of the 
parameter 5 ~ 1 {Vst c>c) corresponds to an isolated 
checkerboard lattice giving the energy gap Egap ~ 0.02. 
However, for infinite staggered potential, s = 1 the en¬ 
ergy gap Egi =0, and validity of results is uncertain be¬ 
cause one cannot restrict calculations to one band only 
when the gap closes. Also, for s close to 1 the spread of 
three states becomes comparable with energy gap, there¬ 
fore their quasi-degeneracy is not visible. For smaller val¬ 
ues of a parameter s, a region with an increased energy 
gap appears (a red color area in Fig. 0, with the largest 
energy gap Egap ~ 0.08 for A ^ 0.1 and the parameter 
s « (0.3, 0.7) {Vst ~ (2.0, 8.0)). Thus, an interaction with 
extra sites, along with opening a single-particle gap, sta¬ 
bilizes the FCI phase. Interestingly, the maximum values 
of many-particle gap coincide with the white line in Fig. 
[2fb) - the transition between indirect and direct gap. 

In Fig. [3] on the right a phase diagram for a phase 
(p = 7r/4 is shown. There are no significant qualitative 
differences comparing to results for p = 7r/2. Quanti¬ 
tatively, the magnitude of many-particle gap is smaller 
than for p = 7r/2. Also, the region of increased gap is 
slightly bigger than for p = 7r/2, because the topological 
phase transition occurs at Vst = V^A instead of Vst = 2A. 

In Fig. m a phase diagram for p = 7r/4 with flat- 
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FIG. 3: A map of the energy gap between 3-fold ground 
state degeneracy and the fourth state for non-flattened 
middle band with (a) phase p = 7r/2 and (b) phase 
(j) = tt/4 (right) as a function of a parameter A and 
a staggered sublattice potential Vst parametrized by 
Vst = Atan(s7r/2). Interaction strengths are Vnn = 
1.5, Vat AT AT = 1. The white line denotes the single-particle 
topological phase transition for Vst = 2A for 0 = 7r/2 and 
Vst = a/ 2A for 0 = 7r/4. 


tened middle band is shown. This corresponds to a map 
of the single particle energy gap Egi from Fig. [2fd). 
Within a major part of a parameters range, the energy 
gap is approximately constant and larger in comparison 
to the energy gaps for non-flattened bands from Fig. [H 
with maximum of Egap ~ 0.085. The single-particle gap 
Egi remains open in the limit Vst oc. Finite value of 
many-particle gap in this limit agrees with earlier results 
for the checkerboard model^*^. No gap closing for finite 
Vst shows that FCI on the Lieb lattice with additional 
hoppings is adiabatically connected to that on checker¬ 
board lattice. A decrease of the energy gap Egap is only 
seen for A ~ 0.1 and close to a single-particle topological 
phase transition {Vst = A) marked by a white line. 

In fig [5] we show the dependence of the energy gap 
between 3-fold ground state degeneracy and the fourth 
states on interaction parameters for fixed A = 0.2 and 
Vst = 2. In general, the energy gap scales approximately 
linearly with an interaction between next-nearest neigh¬ 
bors Vnnn (an interaction between particles occupying 
blue color sites in Fig. [6l^a)) and only slightly depends 
on Vnn (an interaction between particles occupying sites 
with different colors in Fig. [6ja)). This is related to the 
fact that for this choice of parameters the states from 
the middle band are in 98% localized within the sublat¬ 
tice forming a checkerboard lattice, blue color sites in 
Fig. El^a). 


gap 
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FIG. 4: A map of the energy gap between 3-fold ground 
state degeneracy and the fourth state for a middle band 
flattened using additional hoppings t 2 and ts, for a 
phase (j) = 7r/4, as a function of a parameter A and 
a staggered sublattice potential Vst parametrized by 
Vst = Atan(57r/2). Interaction strengths are Vnn = 1-5, 
Vnnn = 1- The white line denotes the single-particle 
topological phase transition for Vst = A. 



FIG. 5: A map of the energy gap as a function of inter¬ 
action parameters Vnn and Vnnn for A = 0.2, Vst = 2. 


B. Identification of FCI phase 

We would like to confirm whether red regions with 
larger energy gaps in Figs. [3] and |4] correspond to FCI 
phase. Thus, for chosen parameters from this region, 
A = 0.2 and lAt = 2, we investigate signatures of 1/3 
Laughlin-like state. Fig. [6fa) shows momentum-resolved 
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FIG. 6: a) Momentum-resolved low energy spectra for 
systems with different sizes given by x Ny) for pa¬ 
rameters A = 0.2, Vst = 2, Vnn = 1-5 and Vnn = 1- The 
energy is rescaled so that ground state energy is set to 0. 
The momenta of 3 quasi-degenerate states agree with a 
counting rule for FCI. (b) Spectral flow upon flux inser¬ 
tion for 4x6 lattice. The 3-fold degenerate ground-states 
flow into each other and do not cross with higher energy 
states. The inset shows magnified view of the ground 
state manifold evolution. 


energy spectrum for different torus sizes. The energy 
spectra are plotted with respect to the ground state en¬ 
ergy at = 0. We find that for each system size we have 
3-fold quasi-degenerate ground state, whose momentum 
counting corresponds to that obtained from generalized 
Pauli principle^. In the case of x Ny = (4x6) this 
correspond to total momenta of three quasi degenerate 
ground states for momenta (Kx^Ky): (0, 0), (0,4), (0,8). 
The electron density of the ground state manifold is al¬ 
most uniformly distributed within sublattice B, as ex¬ 
pected for the incompressible liquid. Small variations 
can be attributed to finite size effects. However, due 
to localization of single-particle wavefunctions on sublat¬ 
tice B, the sublattice A is significantly less filled. In Fig. 
[6](b) the spectral flow upon magnetic flux insertion for 
(4 X 6) torus is shown. The 3-fold degenerate ground 
states do not intersect with higher states. Three states 
flow into each other and return to themselves after inser¬ 
tion of three magnetic fluxes. This no-mixing property 
of the ground state manifold with higher energy states 
is necessary but not sufficient to prove the existence of 
a Laughlin-like phase. Thus, we analyze quasihole ex¬ 
citations from this state^i^. Figure [7| shows quasihole 
spectra for = 7 electrons on (4 x 6) torus (3 quasi¬ 
holes). In this case, 12 quasihole states per momentum 
sector for Laughlin-like excitations is predicted. This is 
indeed observed in Fig. [71 Similarily, the results for 5x5 
torus filled by 8 electrons (one quasihole) also obeys the 
counting rules. The spectrum is divided into two parts 
separated by a clear energy gap, with 12 quasihole states 
per momentum sector below the gap. Thus, our results 
strongly suggest the presence of FCI in this system. 
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FIG. 7: Momentum-resolved energy spectrum for N = 7 
electrons on (4 x 6) torus for parameters A = 0.2, 
Vst = ‘^yNN = 1-5 and Vnn = 1- The number of states 
below the gap starting around E = 0.25 is 12 for each 
momentum sector. This is in agreement with counting 
for Laughlin quasihole states. 


IV. CONCLUSIONS 

In summary, we have analyzed the transition between 
a checkerboard lattice and a Lieb lattice in the context 
of FCI phase for 1/3 filling of a topologically nontrivial 
energy band. Results were presented for two different 
complex phases, and a model with a flattened topologi¬ 
cally nontrivial band. For the non-flattened bands, the 
additional sites open the single-particle energy gap and 
allow FCI to exist. For a flattened band, they increase the 
single-particle energy gap and stabilize the FCI. The exis¬ 
tence of FCI is proven by topological degeneracy, spectral 
flow and momentum counting, both for exact 1/3 filling 
and systems with quasiholes. The topologically nontriv¬ 
ial character of FCI phase is also seen by the fact that it 
exists only in parameter region corresponding to single¬ 
particle topologically nontrivial band. 
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